Chapter 2 
Problem 2.1 
A 90 N/mm 60 N/mm 
A 
u 

Fal 
A 
) 40 N/mm 


Figure S2.1 


Referring to Figure S2.1 the springs with stiff- 
ness 60 N/mm and 90 N/mm are placed in series 
and have an effective stiffness given by 


1 


k a a ee 
‘1/60 + 1/90 — 


36 N/mm 
This combination is now placed in parallel with the 


spring of stiffness 40 N/mm giving a final effective 
stiffness of 


kot = kı +40 = 76 N/mm 


Problem 2.2 


dxdy 


GT » 


Figure $2.2 


The infinitesimal area shown in Figure $2.2(a) 
is equal to rdOdr. When the circular disc moves 
in the x direction with acceleration ti, the inertia 
force on the infinitesimal are is yrd@drii,, where y 
ids the mass per unit area. The resultant inertia 
force on the disc acting in the negative x direction 
is given by 


d 2r 
I; = f / yiizrdodr = yn Ritz = Müs 
0 0 


where M is the total mass of the disc. The resultant 
moment of the inertia forces about the centre of the 
disc, which is also the centre of mass is given by 


"i sin ddédr = 0 


https: Popa 


In a similar manner we get 
Iy = Müy 


For an angular acceleration Ë about the center 
of mass the inertia force on the infinitesimal ele- 
ment is directed along the tangent and is yr?6d0dr. 
The x component of this force is yr20d6dr sin 0. 
It is easily seen that the resultant of all x direc- 
tion forces is zero. In a similar manner the resul- 
tant y direction force is zero. However, a clockwise 
moment about the center of the disc exists and is 
given by 


27 2. 2 
Mg = [ f yör?dðdr = yn Re a =i 


The elliptical plate shown in Figure $2.2(c) is 
divided into the infinitesimal elements as shown. 
The mass of an element is ydxdy and the inertia 
force acting on it when the disc undergoes trans- 
lation in the x direction with acceleration ti, is 
ytizdedy. The resultant inertia force in the neg- 
ative x direction is given by 


a/2 b/24/1— J1—4a? /a? 
L = | it Vuydydx 
a/2 b/24/1— sf 1-42? /a? 
a/2 
sdh a bo/1 — 4 Pde 


a/2 


b 
= = = Mü, 


The moment of the x direction inertia force on an 
element is yiizydvdy. The resultant moment ob- 
tained over the area is zero. The inertia force pro- 
duced by an acceleration in the y direction is ob- 
tained in a similar manner and is Mw, directed in 
the negative y direction. 

An angular acceleration 6 produces a clockwise 
moment equal to yr?0dxdy = y (x? +y?) Odzxdy. 
Integration over the area yields the resultant mo- 
ment, which is clockwise 


a/2 b/2 of 1-42? /a? of 1-42? /a? 7 
Ig = / Von eae y0 (z? + y’) dydx 


a/2 b/24/1— 1-40? /a? 
ag UTES yO ts 
4 16 16 


The x and y direction inertia forces produced on 
the infinitesimal element are — sin Odxdy and 
7@ cos Odxdy, respectively. When summed over the 
area the net forces produced by these are easily 
shown to be zero. 
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Problem 2.3 


Figure $2.3 


As the center of the pulley moves a distance u, 
the spring is stretched by 2u and the pulley rotates 
through u/R. Referring to the free body diagrams 
shown in Figure $2.3, the conditions of equilibrium 
are 


Quk x 2R+miR+ hd -TR=0 
T+ Mi=0 


Substituting Jo = mR?/2 and eliminating T gives 


(a+) i+abu=0 


Problem 2.4 


3/0" 
Light, flexible 


Light, flexible 3 


1'0" 


Wheel 


u 


Figure S2.4 


1 iy. a 
b= ay x1x (3) ~ 768 


a p-2 
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The displacement at location A, shown in Figure 
$2.4, due to a unit load is given by 


OG g 

— 3E,  3EIL, 

_ 12° x 768 36° x 64 

~ 3x 380x108 § 3x 30x 106 xr 
12? x 36 x 32 
12 x 106 xr 

= 0.01475 + 0.01056 + 0.0044 

= 0.0297 in. 


RL; 


Aa mereka 


Equivalent spring stiffness = 


.— 10 _ Ib.s? 
Mass = 3864 — 0.02588 aa 


The equation of motion is given by 


1 = lb. 
0.0297 ~~ 33.67 in. 


0.02588u + 33.67u = 0 


or 
ü + 1300u = 0 


Problem 2.5 


ML, + R)Ö 


cd 


kL,0 


Figure $2.5 


Taking moments about B shown in Figure 
$2.5, and noting that Io = MR?/2 


MR?.. ; 
M (Lı +R) 64 í 6 + cd’ + kL20 = 0 
or 
3 : ; 
M (124 3R? + 211R) Ü+ ed + k30 =0 


Problem 2.6 


kL@ 


Figure $2.6 


Taking moments about A shown in Figure 52.6 
L? m ; 
(“= + us?) Ö + cd?6+kL70 =0 


When the flywheel is braked, an additional rota- 
tional inertia term will exist. Moment equilibrium 


gives 
2 2 .. É 
= +M G Fa +) 6+ cd?Å + kL?0=0 
Problem 2.7 
Boundary conditions 
sal 
dy r=0 
2r i 
dx 
p= l 
dy qL 
EI—~ =0 
dx? 
where 
dp dp 1 1 
ane a ems (86° — 15? + 6€) 
d? 1 
a Sr (24€? — 30€ + 6) 


All four boundary conditions are satisfied 
Ü 
m* =| m {w(x)}° dz 
0 
1 
= mu f {w(€)}° dé = 0.03016mL 
0 
L 2 2 
d(x) 
= EI 
k f { da \ da 


1 2 
_EI [TE IE E 
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Problem 2.8 


Input 


Figure $2.8 


= AlE) = Aı + (A2 — A1) E 
dy 1 en 
LC727 ME) = PAL) 


where p is the mass density. 


L 
m* = f wie) {p(x}? de 


II 


1 
Í iA Cy Ae 
0 
(4£2 + e4 — 46%) de 
= a (5A; + 11A2) 


wf EA(z) = i dz 


0 dx 


1 
-7f {4 + (A2 — A1) E} 
(4 + 4€? — 8€) dé 


E 1 
= 7 (4+34) 


L 
p =- f M(x) tigp(a)dx 


* 


= -f pLitig {Ay + (A2 — Ar) £} (26 — £”) dé 


0 


ane gall 5 
pLiig (GA + 34) 


Problem 2.9 


Refer to Figure 52.9 where the inertia force 


on an infinitesimal element is shown as moe dx. 
The total displacement u% = u + ug, where u is 
the displacement relative to the ground and ug is 
the nacen of E Po The inertia force 
NnaMmics—o S 


uctures-solution/ 


thus becomes mM (Z(t)w(x) + üg) dx. The first term 
within the parentheses produces the generalized 
mass m* while the second term produces the gen- 
eralized force p* as follows 


r aon fs “e 


The generalized stiffness is given i 


EO 


The geometric stiffness produced by the axial force 
S(x) = (L — x)Mg is obtained from 


kg = [ S(x) [w'(a)] ae = sing 


The equation of motion including the gravity effect 


is 
EI 
ae + C = smg) z= < ai 


T40 L ë 8 E 
x 
A 
B(x 
pa 4 
mow dx A dx 
Ê w ut j 
| ae 4 W(X) + Ug 
i 2 
i P(X) M= pie = El z(t) p(x) 
i x ot 
S(x) =(L—x) mg 
v 
Ug 
Figure S2.9 


Problem 2.10 


gaa TE Pyp T’? TE 
55i d2 PTI 


L _ 
5 ~ hey NE mL 
m* = msin T255 


EIZ 
ps 253 


p* -f Fô (a — vt) sin Z dz 
0 L 


Tut 
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Equation of motion: 


mL wtEI . nut 
z = F sin — 
2 213 L 
F 
Sy 


Figure $2.10 


Problem 2.11 


The mass matrix M and the stiffness matrix 
K are given by 


20 0 
M=]0 2 0f kips?/in. 
001 


1000 —500 0 
—500 750 —250 
0 —250 250 


K= kips/in. 


With yT = [1 2 3], the generalized mass 
m* and generalized stiffness k* are obtained from 


m* = TMy = 19 kip s?/in. 
k* = 47 Ky = 1250 kips/in. 
The generalized force is obtained from 


1 


—U'M | 1 | ü = —1170sin(6zt) kips 


The equation of motion is given by 


192 + 12502 = —1170 sin(6rt) 


Problem 2.12 


Taking moments about A shown in Figure 
$2.12 


ml?6 + 2ka?0 + mgl sin @ = 0 


For small vibrations sin 0 ~ 0, hence 


ml?6 + (2ka? + mgl) 0=0 


Figure $2.12 
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Chapter 3 Collectively, the four equations can be expressed in 
the matrix form 


Problem 3.1 
M*ü + C*i+ K*u=0 


u” = [q aq. q3 qa] 
q2 4% W 
W, r g 
Wr 
M* = ge 
|}<—— Q—_>}<—_ Q—__ >| g 7 
g 
T 1 o 2 0 —1 -1l 
BELTS aT 0 2 1—1 
w w Pa 
SS a a d 
kr le -1 24. 0 1 
WT: dd 2 0 —1 —1 
0 2 1 —1 
E a 
| E ip 0 
cs(q3—q1+q28) C5(da— 91-928) =i =l 0 i+ = 
| | kolga = q1 + qaa) | | k5(Q4- q1 — 928) 
Problem 3.2 
krq3 | | wW.. krq4 w. 
g” | | g” mg 
A, do 
w'y 
g” h 
Cs(q3—q1+q28) cs(q4— 91-928) 
+k,(q3— q1 + 92a) + ks(q4— 91 — 928) 
Figure $3.1 p i oa hy 
The degrees of freedom are identified in Figure dúo f eh Mas DON elo eo) 
$3.1. The free-body diagrams are also shown. For 
force balance on each of the front and rear axles b b 
w.. , , ` P 
ge + cs (da — 1 — Goa) Figure $3.2 
+ ks (q4 — qı — 2a) + krqa = 0 Coordinates ut and 0 measure the vertical mo- 
w. , , ; tion and rotation of the base. The forces on the 
Täs + cs (ds — di + 42a) structure are identified in Figure $3.2. Assuming 
f that u’ is measured from the position of equilib- 
+ ks (q3 — qı + q2a) + kras = 0 


rium under gravity load, the condition of vertical 


For equilibrium of vertical forces on the vehicle force balance gives 


body 
ğı + cs (43 + da — 241) mii’ + 2cù + 2ku = 0 (1) 
g 
+ ks (qs + q4 — 2m) = 0 or 
Taking moments about the mass center of the ve- mü + 2ct + 2ku = —miig 
Hiele baay Taking moments about A and assuming that 0 is 
We small 


a ğ2 + csa (ġ3 — ġa + 2aġ2) 


7 Oe er eee | 
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Problem 3.3 


z 

II 

(= 
ošo 


Figure S3.3 


To obtain the first column of stiffness matrix, 
give a unit displacement u, with v and w each zero. 
The resisting forces provided by frames A and B 
and the external forces required to maintain the 
displacement are shown in Figure $3.3. For equi- 
librium kay = 2k1, kay = 0, k31 =0. 

For u = 0, v = 1 and 0 = 0 equilibrium gives 


koa 


k12 = 0, k22 = kə, k32 BA 


In a similar manner, for u= 0, v=0, and 0= 1 


b? a? 
no) eee ey eal 
33 7 + k2 i 


On assembly, the following stiffness matrix is ob- 
tained 


2k 0 0 
K = 0 ko Eza 
0 koa krb? re kaa” 
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Problem 3.4 


For frame A 


uy = 1202 
us = 1264 
or 
0, 
uv} |0 12 0 0 A 
ul |0 0 0 12 03 
04 
Hence 
0 12 0 0 
tai 0 0 A 
Ky, = Tik,T a 
0 0 0 0 
0 540 0 —240 
=n 0 0 0 0 
0 —240 0 240 
For frame B 
Kp=K, 
For frame C 
1 —15 0 0 
te 0 0 1 —15 
860 —12900  —360 5400 
K- -— —12900 193500 5400 —81000 
C™ | 360 5400 360 —5400 
5400 —81000 —5400 193500 
Assembly gives 
860  —12900 —360 5400 
K- —12900 349020 5400 —150120 
E —360 5400 360 —5400 
5400 —151020 —5400 262620 
Mass of slab m = 100%30x24 _ 2,236 Kips’ 
Mass of moment of inertia 
a? + b? 302 + 24? 
Ip =m T = 2.236 x — 5 
= 275.03 kip.ft.s? 
Hence 
2.236 0 0 0 
B 0 275.03 0 0 
M= 0 0 2.236 0 
0 0 0 275.03 
_ 30x12 


x posin Qt = 0.360posin Qt, p =0 , 
of-structures-solution/ 


30 x 6 
1000 
The equations of motion are 


p3 = x posin Qt = 0.180p9 sin Qt, p4 =0 


MÖ +K =p 
where 
67 =[0, 02 03 04] 
p’=[p1 po ps pa] 
Problem 3.5 


uy; Up ug 


kyo Koo kao 
Figure $3.5 


For the crane, apply a unit force at u with 
restraining forces at uz and ug as shown in Figure 
$3.5. 

1 480 


2 
A, =2 = x — x4 -x4 
1 « (5x a ee s0) 


= 0.0683 in. 


Force required for unit displacement = 14.65 kips. 
Force for unit displacement of crane girder at uj, is 


48E 
= 426.04 kips 
(240)3 6.04 kips 


Thus 


ki = 426.04 + 14.65 = 440.7 kips/in 
ko, = —2 x 14.65 = —29.3kips/in 
kzı = 14.65 kips/in 
ASE Ty 
22 — (960)8 


= 58.59 kips/in 
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By symmetry 


k33 = 440.7 kips/in 
k23 = —29.3 kips/in 
kı3 = 14.65 kips/in 


Assembly gives the stiffness matrix 


440.7 —29.3 14.7 
K= | -29.3 58.6 —29.3 | kips/in 
14.7 —29.3 440.7 
Also 
13750 kip.s? 
mh eg a 
™33 = 1711 
= (27500 +7500) _ kip.s? 
nam 386.4 SHOE 
0.0356 0 0 
M=1| o 0.0906 0 
0 0 0.0356 


The equations of motion are 
Mii + ku = 0 


where 


Problem 3.6 


Referring to the free body diagram Figure $3.6 
and taking moments about A and B respectively 


. , _ ka? 

mı Lü + cLùi ia T (ui = u2) =0 
. , ka? 

mMəLüz T BA (uz = u1) =0 


The equations of motion are 


Mü + Cù + Ku=0 


where 
My E 0 
=| ma] S= le o 
ka? 1 —1 
k-i i a= [ u] 


Figure $3.6 


Problem 3.7 
4 


m iĝ 


mg mg 
Figure S3.7 

The free-body diagram for the double pendu- 
lum is shown in Figure $3.7 when the two pendu- 
lums are displaced from the vertical by angles 6, 
and 62. The forces include reactions at the support 
hinges A and B, the inertia forces acting on masses 
mı and mg, the gravity forces mıg and mag, and 
the equal and opposite spring forces ka(@2 — 61). 
Taking moments about A and B we get the follow- 
ing two equations of motion 


ml? — ka? (02 — 61) 
mol? + ka? (02 — 61) 
Noting that for small 0, sin 0 = 0, the two equations 
reduce to 
i++ (“) (6, — 02) + 26, =0 


My l l 


myglsin 6; = 0 


mggl sin 62 = 0 


fi, + © (7), ® 61) + 26) =0 | 
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Problem 3.8 


Referring to Figure $3.8(a), the mass matrix is 
given by 
M= k 0 | 


0 m 


(+) 


ly 4 L4 
SL (d) (e) 
Figure $3.8 


The moments produced by unit forces acting 
along d.o.f. 1 and 2 are shown in Figures $3.8 (b) 
and (c), respectively. The flexibility influence coef- 
ficients are given by 


ayy = i E 2 Sie be 
EI 2 EI 3” T 3EI 
MgMg (Lf/2)* L L L 
= | ay ae a m a 
71 
~ 94 EI 
mımə 1 L L L 
m= ap Sna a ae 


The flexibility matrix is thus given by 


_ 144 EI 
~ 5 pB 


To determine the force along say d.o.f 1 that is 
equivalent to the applied force P, we impose dis- 
placements du, and 0, respectively along the 2 d.o.f. 
The forces required to impose these displacements 
are given by 


[dur] _ [42/5] EI 
fi -x| 0 | = Pa re 


The moments produced by these forces are shown 
in the M diagram Figure $3.8(d). To determine 
the deflection under the applied load P, we apply 
a unit load at the point of application of P and in 
the same direction as P. The m diagram produced 
by the unit load is shown in Figure $3.8(e). The 
deflection is then given by 


Mm 
AU, = | ri 
slig 8 LL 24 
“| 2” SBE A a 3 a4 
36 L L L 
SEI 4° 4° 8 
ig Oe ae 
2” 5EI 4 
1 42  L]EI 
a =| From 
1 9 2l 9 
| 160 mt asl dui = g 


Applying the principle of virtual displacement we 
get 


Qıôuı = PAU, 
which gives 
9 
Qı 160 


In a similar manner, to determine the equivalent 
force along d.o.f. 2 we apply the displacements 0 
and dug along the two d.o.f.. The forces required 
to impose these displacements are given by 


_ 4 [ 0 ]_ [36/5] EI 
=K H = Ha Tsou 


Using the previous results the corresponding dis- 
placement under load P is obtained from 


48 15 9 36 (21 


Using the virtual displacement principle we get 


17 


The equations of motion are thus 
7 1 


0 m|] |i 5 3 |, | lu 
4 3 
9- 
160 
ofa 
40 


w(x) =1 — cos S7 
37x 
Wo(x) =1 — cos y 
iz) = T sin 5 
(eo) = st sin E 


TL 


TSŽ 
r) = (=) cos SF 


The mass influence coefficients are given by 


L 
mi = I m (yı (£)} dz = 0.2268mL 
0 
L 
M22 = | M (Ya(£)) dz = 1.9244mL 
0 
L 


mi2 = M21 = ji mw (a) Wo(x))dax = 0.5755mL 
0 


The stiffness influence coefficients are given by 
L 4 
EI 
ky = | EI (by (x))?dx = —— 
n= f BEGG) Pd Fs 
TEI 
328 


L 
k22 = | EI (2"(a))?da = 81 
0 


L 
kız = kay = J EI (w(x) we!" (ax)) da =0 
0 


The components of the effective load vector are 
given by 


i7 
pı = -müs f ypı(x)dz = —0.3634mLii, 
0 


L 
p2 = -müs f| bo(a)dx = —1.2122mLiig 
0 


The influence coefficients in the geometric stiffness 
matrix are obtained from 


L = 
iers A mg(L — 2) (Wi(2))? de = 0.30997 
0 
L = 
ka22 = | mg(L — x) CAN dx = 5.3017 
0 


L = 
ijas | ing(L — 2) (wh(a)v$(a)) de = 0.7572 


Qz = -—P _ 
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The equations of motion are 


0.5755 1.9244 
: mET[1 0]_ [0370 0.75 z2 
3213 |0 81 I| 0.75 5.302 zo 
0.3634 
1.2122 


_ eo e fl 
mL if 
£2 


= mkü, | 


Problem 3.10 


6El, /Lé 


6El, /Lé 


4El,/Lp 2El,/Lp 


6El,/L2 t 


4El,IL, 


S 6El,/L2 J 
2El,/L, 
6Elp/LE 6El, /LE 
A 2El/Lp 6EI,/L2 
+ { > 


4Ely/Lp 4EI,/L 


2EI/L, A 
=e Sr 
6El,/L2 
12El/LÈ 12El/L? 
6El,/L2 
6El,/L2 6El,/L2 
12El,/L3 12El,/L3 


Figure S3.10 


The first column of the stiffness matrix of an 
individual frame is obtained by imposing displace- 
ments u1 = 1, ug = 0, ug = 0. The external forces 
required to maintain the displacement pattern are 
identified in Figure $3.10 and are given by 


kı = 4El,/ Ly +4811. 
ka = 2EIp/ Ly 


= 2 
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The other columns of the stiffness matrix are 
calculated by a similar procedure. The assembled 
stiffness matrix is given by 


4El, | 4El. 2El, 6EI, 
| Ly Le Ly L2 ] 
2EI, 4B, | 4EI, 6EL 

Ly Ee LZ 


k= Ts 2 
6EI. 6EL, 24ET, 
L L2 Ls 


On substituting the values of E, Ia, Ie, Ly and Le 


A B 
120.0 30.0 0.900 
A 
K = 10° 30.0 120.0 0.900 
B 0.9 0.9 0.036 


Degrees of freedom 1 and 2 are condensed out from 
K to give 


K* = (Kgg — KgaKy4Kazs) = 252 kips/in. 


For frame A 


U3 = 0, + ads 
01 
=[{1 0 a] } %> 
03 
= T40 
K4 = T’K*T 
1 0 120 
= 252 0 0 0 
120 0 14400 
For frame B 
Tg=(0 1 a] 
0 0 0 
Kp = 252 | 0 1 120 
0 120 14400 
For frame C 
T.=[1 0 —a] 
1 0 —120 
K. = 252 0 0 0 
—120 0 14400 
Assembly gives 
K=K4+Kg+K. 
2 0 0 
= 252 | 0 1 120 


0 120 , 43200 
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Slab mass m = 19x20x20 = 0,1035 Sips For the deck 
Mass moment of inertia is given by L 5 
m= f mia) de 
2 4b? 2402 + 240? 
gem sii o mL 
= 993.79 f r ie 
k* = EI TOJ dz = 
213 
0.1035 0 0 0 
M=] 0 0.1035 0 p* = py(zi) 
0 0 993.8 


where x; = vt and p = k(uy — N) + c(t — ù) — me — 


The equations of motion are mg. Thus, we have 


MÖ + Kd =0 


pe = k {us = n} a + c{tty — Hh} a — miia -mga 
where gf = [04 A 03 | 
Equation of motion of the deck is given by 


Problem 3.11 and Problem 3.12 


* 


m*ž + k*z=p 


m, _ u On combining with the equation for sprung mass 
and substituting for 7, ù and 7, the complete set of 
—v equations can be written as 
k E e 
; My 0 Uy 
Up = Z Sin(7x/L) | 0 mL m | | 5 
M, i . 
$ [ex ilg —ca thy 
Be ca? are | Zz | 
L k —ka — c08 Uy 
i m 7 2 2,2, nEI 
—ka ka’ + c0aB— m0 a" + Sy z 
m, Üy 
ee My 
Ku — n(vt)} H {uy — 
mg 
E EE T Problem 3.13 
mi 1} mg 2 
1 f 
pa VOU anas 
Figure S3.11 386.4 
10 tbs s? 
m= — 
=a rE nT 3864 i 
x) = sin — up = zsin — 
Do L TET 
Referring to the free body diagram for the sprung GJ, = © x 12 x 10 = 1, 178, 100 


mass shown in Figure $3.11 2 


EI, = — x 30 x 10° = 1, 472, 620 
mutiny + e(tty — Å) + k(uv — 1) = -mvg gai 


where The degrees of freedom are identified in Fig- 
= za ure $3.13. To derive the flexibility matrix, unit 
ù = z200 + ża loads are applied as shown and the corresponding 


ij = —z0?a + 2206 + ža moment and torque diagrams are drawn. Flexibil- 
ity influence coefficients are obtained by integrating 


6 the 
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netp 


Figure $3.13 


1 12 2 
= er a 
1 367 2 12 
a EE x z X 36 Gy x 36 x 12 
= 0.02970 
12 1 x 36 
=1 14 1 
a22 x El; x GJ, x 
= 3.3775 x 107+ 
1 12 
a12 = oa Oa, x 36 x 1 
= 2.2099 x 1073 
_ 4n—4 | 297.00 22.099 
Aa Bee 3.378 
A eS 65.60 —429.2 
eee e- 5768.1 


Equations of motion are given by 
0.02588 0 ü 
0 0.3170 0 


4 65.60 —429.2 uj |0 
5768.1 0j jo 


—429.2 
Problem 3.14 


Considering that the members are axially rigid, 
end C is fixed and end A can only rotate but not 
translate in either X or Y direction, the structure 
has only two degrees of freedom as identified in Fig- 
ure $3.14 . For segment BC the displaced shape 
corresponding to unit rotation at B can be taken as 


waa(-3) 


where x is measured from B along BC and L = 


2. z miş the length segment. The stiff 
19 gi Smeh. om/ produ 
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k22 is given by 


L 
ae J Ely (a) by (a)de 
0 
_ 4EI _ 4x 20000 x 20 
~~ L 2.5 
= 64.00 x 10* N.m. 


6 kN/m 
CVO VEE VEY 
— 

ici? ay om 
3.5m 
A 

25m >| 

Figure $3.14 
Mass m 

Due to self weight 138.0 kg/m 
Due to superimposed load ie 611.6 kg/m 
Total 749.6 kg/m 


Using this value of total distributed mass 


L 
m2 = | my (ax) v1 (a) da 
0 


_ 749.6 x (2.5)8 
7 105 
= 111.5 kg.m? 


For segment AB mass and stiffness matrices are de- 
rived from shape functions 


uzeli- 


where x is measured from A along AB, and L = 
3.808 m. 

Mass variation is given by m(x) = ma (1+ 2) 
where ma = 69 kg. The mass matrix is obtained 


from 
f 


m(z)pi(x)Yz(x)dx 


Mij = 


structures-solution/ 


which gives 


Man = mA 1128 —9L3 
AB 340 |-9L3 131° 
Substituting ma = 69 kg/m and L = 3.808 m 
49.89 —40.82 
Mas = e | 
The variation of cross sectional flexural rigidity is 


given by 
EI(x) = El, (1+ Z) 


The stiffness influence coefficients are obtained 
from kij = i EI (x)y;" (x); (£)dx which gives 


SEI 3EI TEI 
kı = ~“ kı2 = ko = T 22 = A 
Substituting 
EI, 20000 x 10 F 
2 oe SOAR n 
4 [26.26 15.76 
rae be 36.76 | 
Adding the contributions from AB and BC 
_[ 49.89 —40.82 ; 
as e- n ce 
4 [26.26 15.76 
pret ee 100.76 | “™ 
The equations of motion are 
Mu+ Ku=0 
where u? = [u, u2] 


Problem 3.15 


Figures $3.15 (a) through (f) show the dis- 
placed shapes produced by unit displacement ap- 
plied at one of the d.o.f at a time. The figures 
also identify the moments and forces required to 
maintain such displacements. The stiffness matrix 
is assembled one column at a time. Thus, the first 
column notes the forces required along d.o.f. u1, u2, 
01, 02, 03, 04, 05, and 0g when a unit displacement 
is imposed against u1. 


24EI 12EI 6EI 0 6EI 0 
T3 T3 T2 T2 
12E 24EI 0 6EI 0 6EI 
EB 8 P i? 
6EI 4EI 2EI 
= 0 Ea a 0 0 
6EI 2EI 8EI 2EI 
0 i? l l l 0 
6EI 2EI 8EI 2EI 
i? 0 0 l l l 
0 6EI 0 0 2EI 4EI 


12 l l 
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Figure $3.15 


We partition the matrix, such that the first 2 
by 2 partition, comprising the first two rows and 
the first two columns, is related to d.o.f. uz and u2, 
while the last 4 by 4 partition, comprising the last 
four rows and last four columns, relates to d.o.f. 01 
through 64. We then have 


EI[ 24 -12 
Kaa= | to 24 | 
420 0 
2 ee 3B: ip) 
Kes= T |p 2 8 2 
002 4 


Kga = KI 


ures-solution/ 


Static condensation of the 0 d.o.f. gives the 
following 2 by 2 d.o.f. 


K* =bfK 44 —KapkKgp-1Kpa 


9.6 a | 


= 3 
Bii By 9.6 


The equation of motion is given by 


m 0] [ti], ZTEI J 9.6 84] [u 
0 m| | te T3 |-8.4 9.6 | | us 


https://gioumeh.com/product/dynamics-of-structures-—solution/ 


